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Instability in Shear of Simply Supported Square Plates 

With Reinforced Hole 

By Wilhelmina D. Kroll 

The critical shear stresses were determined for five elastic, simply supported square 
plates with central circular holes of diameter one-eight and one-fourth the length of a side. 
Comparison of the numerical results obtained with those computed for plates without holes 
shows that, although an unreinforced hole may cause a large reduction in the critical shear 
stress of the plate, reinforcement of the hole by a circular doubler plate causes a substantial 
increase in the shear buckling load. 



I. Introduction 

Holes must frequently be cut in the stressed 
skin surfaces of airplane wings or fuselages to 
provide access to the interior of the wing or fuse- 
lage. To prevent weakening of the entire struc- 
ture thereby, the hole is usually reinforced by a 
circular doubler plate riveted to the inner surface 
of the sheet. 

The present report is one of a series of studies 
undertaken for the Bureau of Aeronautics, Navy 
Department, in order to understand more fully 
the effect of reinforcements in strengthening and 
stabilizing the structure around ;i hole. In 
reference [l], 1 a theoretical investigation was 
presented of the effect of various reinforcements 
on the stresses and displacements near a small 
circular hole in a plane sheet under uniform ten- 
sion in all directions. Numerical values of stress 
concentration and of radial displacement at the 
edge of the hole for various doubler plates indi- 
cated that the greatest reduction in stress con- 
centration for a given volume of reinforcing 
material occurs when the outer radius of the 
reinforcement is smallest and that the reduction 
in distortion of the hole depends principally upon 
the volume of the reinforcement. 

Tests were made on plates with reinforced 
circular holes, under uniform tension in one 
direction only, reference [1], to check the theoreti- 
cal analysis, and it was found that, if the reinforce- 



1 Figures in brackets indicate the literature references at the end of this 
paper. 



ments were attached to the sheet with two con- 
centric rows of rivets or bonded to the sheet, the 
theoretical analysis gave a sufficiently accurate 
description of the stresses and displacements in 
the neighborhood of the hole. 

In reference [2], results of compressive tests of 
curved panels with circular holes unreinforced or 
reinforced with doubler plates are given, which 
show that curvature has no appreciable effect on 
the median fiber strain distribution around a hole 
and that the* plane stress theory can be used to 
reliably predict the stress distribution in the un- 
reinforced regions of curved panels with holes. 

A theoretical analysis for the stress distribution 
in a flat plate near a reinforced circular hole 
loaded by a pin was presented in reference [3]. 
Results obtained in the theoretical analysis and 
those from a test of a plate of sandwich con- 
struction showed good agreement. 

A method is presented in reference [4] for deter- 
mining the instability under edge compression of 
a rectangular plate with a reinforced circular hole. 
The values of critical compressive stresses of 
square plates obtained by this method indicated 
that, although the buckling stress of a plate is re- 
duced only a small amount by the presence of an 
ii n reinforced hole, reinforcement of the hole causes 
a substantial increase in the buckling stress. 
Reinforcement of a hole, therefore, can be expected 
to stabilize as well as strengthen the structure in 
the neighborhood of a hole when the structure is 
under compressive load. 
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Using the method of reference [4], a study is 
made in this report of the effect of the reinforce- 
ment of a circular hole in a square plate on the 
stability of the plate under shearing loads uni- 
formly distributed along the edges of the plate. 

II. Method of Analysis 

In reference [4], an energy method was presented 
for computing the compressive buckling load of a 
simply supported elastic rectangular plate having 
a central circular hole reinforced by a circular 
doubler plate. This method can also be used to 
determine the critical shearing stress for simply 
supported rectangular plates having a central 
circular hole. 

The integrals I x and I 2 , which must be evaluated 
for the plate are: 
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where h = plate thickness, function of x and y; 
x, y= rectangular coordinates with origin at 
center of plate; 
w=lateral deflection of plate; 
D=Eh?l\2{\-n 2 ), flexural rigidity of plate, 

function of x and y; 
M=Poisson's ratio, 0.3; 
a x = tensile stress in x direction; 
a y = tensile stress in y direction; 
y = shear stress in x and y directions; 
o= shear stress in x and y directions far 
from the hole. 
The stress ratios a x /r m , o-yjr^ and t X2/ /to> just 
prior to buckling may be obtained from reference 
[4] in the following way. If the origin of coordi- 
nates is taken as shown in figure 1, shear stresses 
along the edges of the plate will be equivalent to 
tensile stresses acting at the edge of the plate at 
0=45° and equal compressive stresses acting at 
the edge of the plate at 6= —45°. The constants 
F and K of reference [4] are zero in the case of 
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shear loading. The stresses oy/reo, cre/T m , and 
Tre/rco referred to polar coordinates r, 6 are then 
given as 



a r /r a3 = -(A+3CR*/r*+2DR 2 /r 2 ) 
cos 2 (*+ j) 

<relr m = (A+6Br 2 /R 2 +3CR 4 lr') 
cos 2 (V|) 

t,#/t.-= (A+3Br 2 /R 2 -3CR 4 lr*- 
DR 2 /r 2 ) sin 2 (0+ j) 



(2) 



where r } 6= polar coordinates with origin at 
center of hole (see fig. 1); 

6= angle between radius r to point and 
positive horizontal axis. Positive 
6 measured in counter clockwise 
direction ; 
R= radius of hole; 
A } B, C, Z>= coefficients with different values in 
sheet and reinforced region. 




Figure 1. Subdivision of plate into three portions A, B, 
and C, in evaluating 1 2 

The values of A, B, C, and D may be determined 
as outlined in reference [4]. The stresses in rec- 
tangular coordinates may be computed from the 
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stresses in polar coordinates, eq 2, by use of the 
conversion formulas: 



c x = <7 r cos 2 0+a-0 sin 2 B—Trs sin 26 
<Ty = a r sm 2 0+c0 cos 2 0+T r sin 2d 

T xy = -^ Sm 26+TrO cos 2d 



(3) 



The lateral deflection, w, of a simply supported 
rectangular plate of length a and width b can be 
approximated by the first terms in the trigono- 
metric series 
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After substituting eq 4 into eq 1 and performing 
the indicated integration by a numerical proce- 
dure, the integrals I x and I 2 become quadratic 
expressions in terms of the coefficients a n , a 22 , 
a 13 , .... The critical value of r m , the 
shear stress far from the hole, at which buckling 
of the plate occurs is that value of Tco which re- 
duces to zero the determinant of the coefficients 
of a n , a 22 , a 13 , etc. in the set of simultaneous 
equations 
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III. Numerical Integration 

The evaluation of the integrals in eq 1 over the 
surface of the plate could not be done directly, 
since the plates have a rectangular outer boundary 
and a circular inner boundary and involve a stress 
that is a complicated function when expressed in 
rectangular coordinates. 

The integral I ly eq 1, was taken as 



where 7 Ia = the integral I x for the surface enclosed 
by the outer rectangular boundary 
of the plate, taking the flexural 
rigidity D=D h the value of I) in 
the unreinforced portion of the 
plate ; 

7 16 = the integral I { for the surface enclosed 
by the outer circular boundary 
of the reinforcement taking 
D=D 2 — D U where D 2 is the value 
of D in the reinforcement; 

7i c =the integral I x for tbe surface enclosed 
by the circular boundary of the 
hole taking B=D 2 . 

The double integration necessary to evaluate 
I i a could always be done directly, but this was 
not possible in evaluating J 16 and l u . The inte- 
grals I\ h and I\ c were obtained by integrating 
directly with respect to x and using Gauss' method 
of numerical integration, reference [5], for inte- 
grating with respect to ?/. In all cases except 
case 2, five Gauss points were used in this numer- 
ical integral ion as computations showed that, 
in increasing the number of points used from 
three to five, the value of the critical stress was 
changed only 2 percent. If more than five points 
had been used, the change in the value of the 
critical stress would have been even less. For 
case 2, three Gauss points were used in evaluating 
I Xc because the circular area over which I u was 
evaluated had a much smaller radius for case 2 
than for the other cases. 

The integral I 2 was evaluated as the sum of 
three integrals 



I* = Iu + I»+hc, 



(7) 



where I 2a =ihe integral I 2 for the circular disk 
of the reinforced area (A, fig. 1); 
I 2b =the integral I 2 for the circular disk 
between the outer boundary of 
the reinforcement and the largest 
inscribed circle in the plate (B, 
fie* r 
the integral I 2 for the remainder of 
the plate (C, fig. 1). 
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The integral for each circular portion was deter- 
mined, using Gauss' method of numerical inte- 
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gration, by first integrating numerically in a 
circumferential direction and then in a radial 
direction. The integral for the remainder of the 
plate was obtained also by using Gauss' method, 
first integrating numerically in the ^-direction and 
then in the ^/-direction. 



A typical distribution of Gauss points for com- 
puting I 2 is shown in figure 2 for one-eighth of a 
square plate with a central circular hole. For 
this distribution, the integral of a function F over 
the portion of the plate surface enclosed by heavy 
lines in figure 2 is 



jt 



FdA=j [0.0481444^ + 0.0770310^2 + 0.0481444^3 + 0.0380048^4+0.0571868^5 + 
0.0914988^6+0.0571868^7+0.0727220^8 + 0.1163557^9+0.0727220^0 + 
0.0337158^ + 0.0539452^2+0.0337158^3 + 0.0142967^4+0.0228747^5+ 
0.0142967f 16 + 0.0181805F 17 +0.0290888F 18 +0.0181805F 19 + 0.00842899F 20 + 
0.0134864^2i + 0.00842899F 22 ] 



(8) 
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Figure 2. Distribution of Gauss points for computing 
1 2 for case 3, 

Coordinates of points 



Point 


x/a 


y/a 


Point 


x/a 


Via 


1 


0. 469052 


0. 206523 


12 


0. 257000 


0. 106453 


2 


. 469052 


.321115 


13 


. 213318 


. 178541 


3 


. 469052 


. 435706 


14 


. 234988 


. 020855 


4 


. 384501 


. 352062 


15 


.217954 


. 090280 


5 


. 469978 


. 041709 


16 


. 180908 


.151415 


6 


. 435910 


.180559 


17 


. 186766 


. 016575 


7 


.361818 


. 302831 


18 


. 173228 


. 071753 


8 


. 373532 


.033150 


19 


. 143784 


. 120343 


9 


. 346455 


. 143506 


20 


. 138543 


. 012295 


10 


. 287568 


. 240686 


21 


. 128501 


. 053227 


11 


. 277086 


. 024591 


22 


. 106659 


. 089271 



where F\, F 2 , F 3 , . . . F 22 are values of the function 
i^at the positions shown in Figure 2. 

In each case a sufficient number of Gauss points 
was used to reduce the estimated error to less than 
5 percent. Twenty- two points for each octant of 
the plate were used for cases 1, 3, and 4; twenty- 
eight points for case 2; and twenty-five points for 
case 5. 

As an indication of the adequacy of the numer- 
ical integration methods used, the critical shear 
stress of a square plate without a hole was deter- 
mined by approximating the deflection by 



WX 



iry_ 



w—a n cos — cos - +a 22 sin 
a a 



2*y 



i-KX . 

- sin 
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integrating for I u eq 1, exactly and integrating 
for I 2 , eq 1, both exactly and by numerical 
integration with 22 points in one-eighth of the 
plate as shown in figure 3. The resulting critical 
stresses differed by 0.8 percent. In the case of a 
plate with a reinforced hole, it is probable that 
the more complicated stress distribution and the 
use and prominence of higher order terms in 
the series used for the deflection will cause the 
error to be somewhat higher. 

IV. Convergence of Deflection Function 

The correct value of r^ for buckling of the 
plate is that value of Too which reduces the 
determinant of the coefficients of a n , a 22 , a n , . . . , 
in the infinite set of eq 5 to zero. In order to 
limit the work of computing Too to a finite amount, 
preliminary computations were made to deter- 
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mine which terms in the deflection function, eq 4, 
were most important in evaluating Too. The 
computations were made for the plate of case 
3, table 1. 

Table 1. Dimensions for square plates investigated 



Case 


Length 
of side 



a 

a 
a 
a 


Radius of 
hole, R 


Outer 
radius of 
reinforce- 
ment 

(«) 
(«) 

0.25 a 
.25 a 
. 1875 a 


Thickness 


Reinforced 
region 


Unrein- 

forced 
region 

h 
h 
h 
h 
h 


1 
2 
3 
4 
5 


0. 125 a 
. 0625 a 
.125 a 
. 125 a 
. 125 a 


2.0 h 
1.5 h 
3.4 h 



a Hole in plate not reinforced 
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Figure 3. Distribution of Gauss points for computing 
1 2 for case of plate without a hole. 

Using the distribution of points shown in figure 
4 to carry out the numerical integration for I 2 , 
the critical shear stress r m was 18.78 Eh 2 /a 2 using 
only the first two terms of eq 4, 
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Figure 4. Distribution of Gauss points for evaluating 7 2 in 
study of convergence of deflection function. 

Too was 18.10 Eh 2 /a 2 , a decrease of 3.6 percent. 
Using 

ttx 7rty . . 2ttx . 2iry . 

w=a,u cos — cos — +#•» sin - sm l + 
a a a a 

3tx Wiry 
a 33 cos - cos —z-) 
a a 

Too was 17.46 Ek 2 /a 2 , a decrease of 7.0 percent. 
Using 

TTX 



w=a u cos 



7rt/ , . 2wx . 2-iry , 
cos — - + a 22 sin sm ■ — - + 



a a 

. 2irx . 4iry . . \ttx . 2iry 

a 24 sm — sm — + a 42 sm - sin ■ -> 
M a a " a a 

To, was 17.37 Eh 2 1 a 2 , a decrease of 7.5 percent. 
Using 

ttx iry . . 2tx . 2-n-y , 

w=a n cos — cos — +a 22 sm sm — - + 

a a a a 

3 TTX 57r?/ , 57TX 37T?/ 
#35 COS COS — +#53 COS COS > 

a a a a 

Too was 18.50 Eh 2 /a 2 , a decrease of 1.5 percent. 
Using 



TTX iry 



Itx . Ziry 



w=a„ cos — cos —-r-a 2 2 sm sin — - + 

a a a a 



. 4:TTX . 4:iry 

a A a sm — - sm - — ; 
a a 

Too was 18.58 Eh 2 /a 2 , a decrease of 1.1 percent. 
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The use of additional terms of the trigono- 
metric series with the first two terms probably 
would have shown a proportionately smaller 
decrease in the buckling load. Therefore, it is 
believed that the function 

irx tv . . 2irx . 2iry . 

w=a n cos — cos — + ao 2 sin — sin — - + 
a a a a 

SttX StV . . 2irx . \tsv , 

a 33 cos — cos — - + a 24 sm — sin — - + 



a 



. 4irx . 

a 42 sm sm 

^ a a 



a 

2wy 



approximated the lateral deflection of a square 
plate with simply supported edges under shearing 
loads with sufficient accuracy. This deflection 
function was used in the computations of the 
critical shearing stresses. 

After the critical shearing stresses had been 
computed for the five square plates considered in 
this report, the relative importance of the various 
terms in the deflection function was determined. 
Values of a 22 were 25 to 42 percent of a n ; coeffi- 
cients a 33 , a 24 , and a 42 were 0.4 to 5.0 percent of 
a n . It seems reasonable to assume that higher 
order terms would have been even smaller as 
compared to an and that the most important 
terms were used for the deflection in all cases. 

The number of points used in the numerical 
integration for I 2 has an effect on the accuracy of 
the determination of the critical stress. To 
carry out a complete investigation of convergence, 
as more points are used, was not possible because 
of the amount of work involved. To obtain a 
partial indication of the accuracy to be expected, 
however, the critical stress for case 3, table 1, was 
computed using various deflection functions and 
using both the distribution of points in figure 2 
and that in figure 4 for computing I 2 . The results 
are given in table 2. 

Table 2. Values of critical stress obtained from the spacings 
of Gauss points investigated 



Terms in deflection 
function eq 4 


Critical stress 
using points in 
figure 4 for h 


Critical stress 
using points in 
figure 2 for 72 


Differ- 
ence 


All, 022 

011,022,033 

011,022,024,042 

011,022,035,053 


18.78 Eh?-la> 
17.46 Ehya 2 
17.37 EM/a* 
18.50 Eh'-la* 


20.05 Ehya? 
18.56 Eh*/ a* 
18.21 Ehya* 
19.58 Ehya* 


Percent 
7 
6 
5 
6 



It was decided on the basis of these computa- 
tions that the use of the point distribution in 
figure 2 would be acceptable. 

The over-all accuracy of the computations of 
critical stress are estimated to be such that the 
critical stresses are within 7 percent of the exact 
values for the cases investigated. 

V. Results and Discussion 

The critical shear stress was determined for five 
square plates with reinforced and unreinforced 
holes, figure 5. The dimensions for the plates are 
given in table 1. 

The analysis gave the critical shear stress far 
from the hole T m , corresponding to the stress dis- 
tribution for a plate under uniform shear at an 
infinite distance from the hole. The shear stress 
Too is larger than the average shear stress on the 
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Figure 5. Loading of square plates with simply supported 
(Dimensions in table 1.) 
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center cross section obtained by dividing the load 
on that section by ah. A better value of the 
average buckling stress r cr of the plate was ob- 
tained from 

1 P1 + P2 

2 ah ' 



Tcr = 



(9) 



where P Y = shearing force acting on edge of plate; 
P 2 = shearing force transmitted across center 
section, ]/=0. 

The value of Pi was obtained by numerical integra- 
tion and the value of P 2 by direct integration of 
the stresses r xy given by eq 3. Values of T m , 
n=PJah, r 2 =P 2 /ah, and r cr are given in table 3. 

Table 3. Critical shear stress for plates investigated 



Case 


Too 


r, 


T 2 


T fr 











Eft* 

8.44 — 5- 








1 


6.95^ 

a 2 


0.9, '••'!-' 
0- 


«11 ^ 


6.53 - flT - 


2 


a 2 


8.55 -p- 


Eh 2 

8.29^- 

a 2 


8.42^ 


3 


17.47^ 

a 2 


!7 52^ 


Eh 2 
17.01 ^ 


,7.27 » 
a 2 


4 


.».«£ 


10.46^ 

a 2 


9 78^ 
9-78 fl2 


-f 


5 


2,83 f 


28.86^ 

rt 2 


Eh 2 
27.80 „ 


28.33^ 



It was assumed that the boundaries of the 

square plates were sufficiently far from the hole 
so that only a uniform shear was acting on the 
boundaries. Numerical integration of the stresses 
<r x as given by cq 3 on the boundary of the plate 
showed that there was an average compressive 
stress on half the side of the plate and a tensile 
stress on the other half equal to 3 to 12 percent 
of the shear stress. This accounts for the small 
differences in the values of t x and r 2 . 

The shear buckling stress of a simply supported 
square plate of constant thickness is given on 
page 360 of reference [6] (taking Poisson's ratio 
as 0.3) as 

r cr =8.44 Eh 2 /a 2 . (10) 

Comparison of this value, case 0, table 3, with the 
values of the critical stresses for cases 1 and 2 



shows that the unreinforced holes of diameter one- 
fourth and one-eighth the length of a side reduced 
the buckling stresses by 22.6 and 0.2 percent, re- 
spectively. Comparison of the value for case 
with the values for cases 3 to 5 of table 3 shows 
that reinforcement of the hole increased the shear 
buckling stress of the plate over that for a plate 
without a hole from 20 to 236 percent. These 
percentages are to be taken only as indications of 
the order of magnitude since the accuracy of the 
present determination of critical stresses is esti- 
mated to be 7 percent. 

The effect of thickness of the reinforcement on 
the buckling load is seen by comparing cases 3 
and 4, table 3. In these cases, the hole size and 
area of the reinforcing material is the same, but 
the thickness of the reinforcement of case 4 is only 
one-half the thickness of the reinforcement of 
case 3. The 4 buckling load for case 3 is 71 percent 
higher than that for case 4. 

The effect of shape of reinforcement on the 
critical shear stress is indicated by comparing 
cases 3 and 5, table 3. The volume of material 
reinforcing the holes in these two cases is the 
same, but the reinforcing material is concentrated 
nearer the edge of the hole in case 5 than it is in 
case 3. The buckling stress is 64 percent greater 
for case 5 than for case 3. 

Comparison of the values for case 1 and cases 
3 to 5 shows that reinforcement of the hole raises 
the critical shear stress from 55 to 334 percent, 
depending on the thickness and shape of the 
reinforcement . 

VI. Conclusions 

Using a numerical procedure for evaluating the 
integrals for the energy stored in a plate, the 
critical shear stress can be estimated for plates 
with circular holes and doubler plate reinforcement. 

The critical shearing stress of plates may be 
reduced considerably by the presence of unrein- 
forced holes. A square plate with a hole diameter 
one-eighth the length of a side showed practically 
no reduction in critical shearing stress, but a square 
plate with a hole diameter twice as long showed a 
reduction of 23 percent. 
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Reinforcement of the doubler plate type causes 
marked increases in the critical shear stress. The 
critical shearing stresses were 55 to 334 percent 
higher than that for a plate with the same size 
unreinforced hole. 

The volume and shape of the reinforcement has 
a marked effect on the critical shearing stress. 
The critical stresses for two plates with the same 
size hole and same area of reinforcement, but 
different thickness, were compared. The critical 
stress increased from 10.12 Eh 2 1 a 2 to 17.27 Eh 2 1 a 2 , 
an increase of 71 percent, when the thickness of 
the reinforced region was increased from 1.5 h to 
2.0 h. The critical stresses for two plates with 
the same size hole and same volume of reinforcing 
material were compared. The plate with the 
material concentrated closer to the edge of the 
hole had a 64 percent larger critical shearing stress. 
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